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ABSTRACT 

We present (V, /) photometry of two wide (~ 25 x 25 arcmin^) fields centered on the low surface 
brightness dwarf spheroidal galaxies Draco and Ursa Minor. New estimates of the distance to these 
galaxies are provided ((m - M)o{UMi) = 19.41 ± 0.12 and (m - M)Q{Dra) = 19.84 ± 0.14) and a 
comparative study of their evolved stellar population is presented. We detect for the first time the RGB- 
bump in the Luminosity Function of UMi (V^]^^ = 19.40 ± 0.06) while the feature is not detected in 
Draco. Photometric metallicity distributions are obtained for the two galaxies and an accurate analysis 
to determine the intrinsic metallicity spread is performed by means of artificial stars experiments. The 
adopted method is insensitive to stars more metal poor than [Fe/H] ~ — 2.5 and it rests on the assumpion 
that the age spread in the considered populations is small (i.e. the impact of the actual age spread 
on the colors of the RCB stars is negligible). We find that, while the average metallicity of the two 
galaxies is similar (< [Fe/H] >uMi= —1-8 and < [Fe/H] >Dra= ^1-7) the metallicity distributions 
are significantly different, having different peak values {[Fe/HYjj^f^ — —1.9 and [Fe/ H]^°^ — —1.6) 
and different maximum metallicities. We suggest that such differences may be partly responsible for the 
difference in HB morphology between the two galaxies. The intrinsic metallicity 1 — a spread is Ui = 0.10 
in UMi and ai = 0.13 in Draco. We demonstrate that the inner region of UMi is significantly structured, 
at odds with what expected for a system in dynamical equilibrium. In particular we show that the main 
density peak of UMi is off-centered with respect to the center of symmetry of the whole galaxy and it 
shows a much lower ellipticity with respect to the rest of the galaxy. Moreover, UMi stars are shown 
to be clustered according to two different characteristic clustering scales, as opposite to Draco, which 
instead has a very symmetric and smooth density profile. The possible consequences of this striking 
structural difference on our ideas about galaxy formation are briefly discussed. Combining our distance 
modulus with the more recent estimates of the total luminosity of UMi, we find that the mass to light 
(M/L) ratio of this galaxy may be as low as M/L ~ 7, a factor 5-10 lower than current estimates. 

Subject headings: (galaxies:) Local Group - galaxies: dwarf - galaxies: individual (Draco, Ursa Minor) - 
(cosmology:) dark matter - galaxies: structure - galaxies: fundamental parameters - 
galaxies: stellar content 



1. INTRODUCTION 

The Draco (aaooo = 17''20'"19", 62000 = 57°54.8') and 
Ursa Minor (aaooo = 15''09™1P, ^2000 = 67°12.9') dwarf 
spheroidal (dSph) galaxies are the faintest known members 
of the Local Group of galaxies and they are among the low- 
est surface brightness members of the group (Mateo 1998). 
They appear to be dominated by very old (age > 8 — 10 
Gyr) and metal deficient {[Fe/H] ~ —2) stellar popula- 
tions (see Mateo 1998; Carney & Seitzer 1986; Aparicio, 
Carrera & Martinez-Dclgado 2001; Grillmair et al. 1998; 
MigheU & Burke 1999; Dolphin 2002), thus they may rep- 
resent a very early and elementary stage of the evolution 
of the building blocks that may have had a primary role in 
the assembly of the Milky Way galaxy (see Bellazzini, Fer- 



raro & Pancino 2001b). Furthermore, they are reported to 
have the highest Mass to Light (M/L) ratio of any other 
known galaxy (e.g., up to M/L — 300 — 1000 for Draco, 
and M/L - 50 - 100, for UMi, see Kleyna et al. 2001; 
Armandroff, Olszewski & Pryor 1995). Hence their stellar 
content may just represent the handful of baryons trapped 
in a system whose true nature is that of a huge dark halo. 
All these exceptional properties have made these galaxies 
a classical case of study. Nevertheless, serious observa- 
tional problems (e.g., the very low surface brightness that 
requires very wide field photometry to obtain statistically 
significant samples of members stars) have hampered our 
knowledge of these intriguing stellar systems. 

Draco and UMi are twin galaxies under many aspects: 
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they have a similar distance from the eenter of the Galaxy, 
similar masses and luminosities, similar metal content and 
are both devoid of gas (see Table 3 for a summary of their 
physical parameters). In this context, a comparative study 
performed with strictly homogeneous observational mate- 
rial and with the same data reduction / data analysis tech- 
niques may reveal interesting features. Here we present the 
results of a comparative analysis, performed by obtaining 
well calibrated (V,I) photometry of the evolved stars of 
Umi and Draco over a field of view ^ 25 x 25 arcmin^, un- 
der strictly homogeneous conditions (see Bellazzini et al. 
2001, for a discussion of the possible problems associated 
with the comparison of non- homogeneous photometries). 

The plan of the paper is as follows: in §2 we describe the 
observational material, the data reduction process and the 
photometric calibration; in §3 we present the Color Mag- 
nitude Diagrams (CMD) and the results of artificial stars 
experiments; §4 is devoted to the estimate of the distance 
to Draco and Umi; in §5 we study the properties of the 
Red Giant Branch of the two galaxies. In §6 we compare 
the structures of Dra and UMi and we demonstrate that 
the inner region of UMi is significantly structured and that 
its stars are clustered according to two different character- 
istic scale-lenghts. §7 is dedicated to the discussion of our 
results. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Observations 

The data were obtained at the 3.52 m Italian tele- 
scope TNG (Tclcscopio Nazionalc Galileo - Roque de los 
Muchachos, La Palma. Canary Islands, Spain), using Do- 
LoRes, a focal reducer imager/spectrograph equipped with 
a 2048 X 2048 pixels thinned and back-illuminated Loral 
CCD array (gain= 0.97e- /ADU, read-out noise 9.0 ADU 
rms). The pixel scale is 0.275 arcsec/px, thus the total 
field of view of the camera is 9.4 x 9.4 arcniin^. The 
observations were carried out during three nights (March 
19, 20 and 21, 2001), under average seeing conditions 
(FWHM ~ 1.0 - 1.4 arcsec). The first and third nights of 
the nm were photometric. The scientific exposures taken 
during the second night were calibrated indirectly by using 
the data acquired in the photometric nights. 

Each galaxy was sampled with a square mosaic of nine 
partially overlapping fields, covering a total field of fiew of 
25 X 25 arcmiv? . For each field two V and two I expo- 
sures {texp = 120 s in each filter for Dra, and t^xp = 90 s 
for UMi) were secured. 

2.2. Data Analysis 

All the raw images were corrected for bias and flat 
field, and the overscan region was trimmed using standard 
IRAF'^ procedures. Each pair of images, in each band, has 
been averaged, so that the final analysis was performed on 
averaged V and I images. 

The PSF-fitting procedure was performed indepen- 
dently on each V and I average image, using a version 
of DoPhot (Schechter, Mateo, & Saha 1993) modified by 
P. Montegriffo at the Bologna Observatory to read images 
in double precision format. The frames were searched for 

IRAF is distributed by the National Optical Astronomy Observal 
Astronomy, Inc., under cooperative agreement with the National Sci 



sources adopting a b-a threshold, and the spatial vara- 
tions of the PSF were modeled with a quadratic polyno- 
mial. A final catalogue listing the instrumental V,I mag- 
nitudes for all the stars in each field has been obtained by 
cross-correlating the V and I catalogues. Only the sources 
classified as stars by the code have been retained. The 
spurious sources erroneously fitted by DoPhot (as cosmic 
rays, bright background galaxies etc.) have been removed 
by hand from the catalogues. 

Nine different catalogs (one for each sub-field of the mo- 
saic) were obtained for each galaxy. These catalogue were 
reported to a homogeneous photometric system using the 
large sets of stars in common among the various fields. 
The magnitudes of the stars in common among adjacent 
fields were averaged. A homogeneous total catalogue of in- 
strumental magnitudes and positions was finally obtained 
obtained for each mosaic field. 

2.3. Photometric calibration 

The absolute calibration has been obtained from several 

repeated observations of Landolt (1992) standard fields, 
including all the stars listed in the extended catalogue of 
calibrators provided by Stetson (2000). The coefficents of 
atmospheric extintion (Cext) were directly obtained by re- 
peated observations of the same standard field at different 
airmasses. Accurate estimates of the aperture correction 
were obtained with a large number of bright and isolated 
stars. The final calibration equations and Cext are shown 
in Fig. 1. 

The accuracy of the photometric calibration was checked 
by comparison with the independent (V,I) photometry of 
the Draco dSph (kindly provided by Prof. P. Stetson, pri- 
vate communication) , and with the V photometry of UMi 
from Stetson (2000). The results of this test are shown in 
Fig 2 where the difference between our photometry (sub- 
script t.w., e.g. this work) and the one by Stetson (sub- 
script SOO) is plotted as a function of magnitude for both 
galaxies. Unfortunately no previous I band photometry 
of UMi is available. The comparison suggests that the 
global imcertainty in the calibration is lower than ±0.02 
mag in each passband, over the whole range of {V — I) 
color sampled. The equations shown in Fig. 1 were ap- 
plied to the catalogs described in §2.2 to produce the final 
calibrated catalogs that are reported in Table 1 and Ta- 
ble 2. In each table we report: (column 1) the identifica- 
tion number, (col. 2 and 3) the V magnitude and its error 
as estimated by DoPhot (see §3.1, below), (col. 4 and 5) 
I magnitude and error, (col. 6 and 7) the position in the 
mosaic in pixels, (col. 8 and 9) the position in the sky 
in equatorial coordinates at the equinox 2000.0, (col. 10) 
the classification of the identified variable stars, if avail- 
able and (col. 11) the name of the identified variable stars 
according to the nomenclature indicated in the table note. 
The photometric errors as estimated from artificial stars 
experiments will be shown and discussed in §3.1. 

3. THE COLOR MAGNITUDE DIAGRAM 

The final (V, V-I) Color Magnitude Diagrams (CMD) 
of the two galaxies are shown in Fig. 3. In both cases 
our photometry reaches the base of the Red Giant Branch 

', which is operated by the Association of Universities for Research in 
;e Pundation. 
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Fig. 1. — The difference between calibrated (upper case letters, e.g., V,I) and instrumental (lower case letters, e.g., v,i) magnitudes for the 
Landolt standard stars observed on March 19 (triangles) and March 21 (circles). The calibration relations (upper panel: V; lower panel: I) are 
plotted (solid lines) and the corresopnding equations are reported along with the root mean square error of the linear fits and the measuredt 
extinction coefficients. 
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Fic;. 2. — DifFcrcncc between independently calibrated photometry from this work (t.w) and by Stetson (2000), for Draco (upper panel: V, 
middle panel: I) and UMi (lower panel: V). 
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(RGB) sXV ^ 22.5. 

The RGBs are well defined and steep, typical of metal 
poor stellar systems. The RGB of UMi appears to be nar- 
rower than the one of Draco, suggesting a smaller metal- 
licity spread or a similar mctallicity spread but a lower 
mean metal content. In the CMD of UMi a hint of the 
Asymptotic Giant Branch (AGB) bump (sec Fcrraro ct al. 
1999) can be noted (at V ~ 19,F - / ~ 0.9) while the 
corresponding feature of Draco is probably hidden by the 
larger extent of contamination by foreground and back- 
ground stars. Note the sharp color cut-off of the distribu- 
tion of field stars at V— / ~ 0.8, corresponding to the Main 
Sequence Turn Off (MSTO) of the Galactic halo/thick disk 
stars (Morrison et al. 2000). Both Eskridge & Schweitzer 
(2001) and Shctronc, Cote & Stetson (2001) foimd that 
some stars lying to the blue of the upper RGB are mem- 
bers of the respective galaxy, either in Draco or in UMi. 
The same authors have also demonstrated that all of these 
stars are not first-ascent red giants, but carbon stars in- 
stead. 

Most of the Horizontal Branch (HB) stars lie in the 
range 19.5 < V < 20.5 in UMi and 20.0 < F < 21.0 
in Draco. Wc have counter-identified 65 variable stars in 
UMi from the catalogue of Nemec, Wehlau & Mendes de 
Oliveira (1988, hereafter N88), that cover our whole field of 
view. In Fig. 3 we have marked the type ab RR Lyrae with 
open squares, the type c RR Lyrae with open triangles and 
the Anomalous Ccpheids with open circles, after the clas- 
sification by the same authors. In Draco 56 variable stars 
were counter-identified in the central ~ 10 x 10 arcmin? 
field, from the study by Baade & Swope (1961). Since an 
explicit classification is not provided by Baade & Swope 
(1961) wc marked all the RR Lyrae variables with stars, 
while the only Anomalous Cepheid identified is plotted as 
an open circle, as in the previous case. However, it has to 
be recalled that, according to Baade & Swope (1961), the 
large majority of RR Lyrae in Draco arc of type ab (see 
also Nemec 1985). 

As already noted by other authors the morphology of 
the two HBs is very diff'erent. UMi has only an handful 
of HB stars to the red of the instability strip, a significant 
number of RR Lyrae variables of both Bailey's types and a 
well populated blue tail reaching V ^ 20.5 and V—I ~ 0.0. 
On the other hand the HB of Draco is well populated in 
its red part and has a sparse blue tail, reaching V ~ 21.0. 
As can be seen from Fig. 3 there is a large number of stars 
in the red part of the instability strip (see also Fig. 4, be- 
low), suggesting that the census of RR Lyrae variables in 
this galaxy is far from complete. A modern CCD search 
for variables over a wide field of view is urged. 

Finally at 21 < < 22.5 and 0.0 < V" - / < 0.6 the 
upper part of a Blue Plume is evident in both CMDs. This 
is a well known feature of the CMD of these galaxies and 
has been preferentially interpreted as a sequence of gen- 
uine Blue Straggler stars, e.g., the result of the evolution 
of binary stars (Carney & Seitzer 1986; Grillmair et al. 
1998; Mighell & Burke 1999; Carreraet al. 2002; Aparicio, 
Carrera & Martinez-Delgado 2001). 

3.1. Artificial stars experiments 

We have performed extensive artificial stars experiments 
in the central field of each mosaic, i.e. the fields with the 



highest stellar density. We take the results from these 
fields as representative of the whole mosaics. More than 
5000 stars per field have been extracted from a Luminosity 
Function (LF) similar to the observed one, and with col- 
ors lying on the observed RGB ridge line of the galaxies, 
following the methods described in detail in BcUazzini et 
al. (2002a, b). The artificial stars have been simultane- 
ously added to the V and I frames, ^ 400 at a time, with a 
minimum distance among them of ~ 40 px, to avoid unde- 
sired interference among artificial stars (see Bellazzini et 
al. 2002b, for details and references). The stars are sim- 
ulated with the observed Point Spread Function and the 
complete data-reduction process has been repeated for all 
the frames "enriched" with the artificial stars. 

In the present context we performed the artificial stars 
experiment to obtain: 

1. A realistic estimate of the photometric error as a 
function of magnitude, by comparing the a priori 
known input magnitudes with the measured (out- 
put) magnitudes. This result is shown in Fig. 4. 
The average error in lower than 0.05 mag in both 
passbands for V < 20.0 and < 0.1 for fainter mag- 
nitudes. 

2. The effect of the "observation plus data reduction" 
process on the measured magnitudes and colors of 
stars lying along the RGB of a Simple Stellar Popu- 
lation (SSP; e.g. a population of stars all having the 
same age and chemical composition; see Renzini & 
Fusi Pecci 1988). This result (shown in the upper 
panels of Fig. 5) will be very useful to discriminate 
the color spread due to the intrinsic metallicity dis- 
persion from the observational scatter. 

3. The completeness factor (C/ = Nrec/Ninpi i-e. the 
ratio between the artificial stars correctely recovered 
and measured and the total number of simulated 
stars) as a function of magnitude. The lower panels 
of Fig. 5 show that, in both mosaics, the complete- 
ness is very high (C/ > 90 %) down to 7 20.5. 

4. THE DISTANCE TO UMI AND DRA 

4.1. Relative distance 

As a first step wc determine the relative distance mod- 
ulus of the two galaxies. Since they two are affected by 
the same amount of interstellar extinction (E(B-V)= 0.03 
Mateo 1998; Schlegel, Finkbeiner & Davis 1998) and they 
share a similar metal content (Mateo 1998), such estimate 
can be obtained by finding the best match between the 
HBs. The upper panel of Fig. 6 shows a zoomed view of 
the HB of UMi. The RR Lyrae are plotted as open stars. 
The stars not recognized as variables by N88 but lying 
in the range of color and magnitude covered by the RR 
Lyrae sampled at random phase are plotted as dots, while 
the stars to the red and to the blue of this region are plot- 
ted as solid circles. The solid line in the instability strip is 
the mean RR Lyrae level (< Vrr >= 19.86 ±0.07) accord- 
ing to N88. The solid heavy line is the ridge line tracking 
the uician of the HB distribution of UMi. The dashed lines 
are displaced by AV = ±0.1 mag from the ridge line and 
roughly bounds all the non- variable HB stars of UMi.. 
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Fig. 3. — Color Magnitude Diagram of 4439 stars in UMi (left panel), and of 4585 stars in Draco (right panel). The Anomalous Cepheids 
are plotted as open circles in both panels. In the left panel (UMi) the ab type RR Lyrae are plotted as open squares and the c type as 
triangles. In the right panel (Dra) all the RR Ly stars are plotted as stars, independently of their type. 



The Draco and UMi dwarf spheroidals 



7 



0.4 

0.2 



-0.2 
-0.4 F 



T 1 1 r 



T 1 1 1 r 



UMi 



_l I I I I I I L. 



T 1 1 r 



_l I L. 



17 



18 



19 



out 



20 



21 



0.4 



3 0.2 

O 



I 



r 



>" -0.2 r 
-0.4 r 




J I I I I I I I I I I I I I I I I I L. 



J I I L. 



18 



19 



20 



21 



22 



0.4 


~ Dra 


1 1 1 1 1 1 1 1 1 1 1 - 


0.2 

















0.2 






0.4 


1 1 1 1 1 1 1 1 1 1 1 1 


I 1 1 1 1 1 1 1 1 1 - 



17 18 19 20 21 



U.4 _ 
g 0.2 ^ 

o; ^ ^ ^ — ^ ^ 

>" -0.2 '- 

18 19 20 21 22 

V 



Fig. 4. — Magnitude residuals from artificial stars experiment! 




Fig. 5. — Upper panels: the input and output CMDs of the artificial stars for UMi (left pair of panels) and Dra (right pair of panels). 
These plots show the effects of the observation + data reduction process on the RGB stars of a Simple Stellar Population. Lower panels: The 
completeness factor (Cy) as a function of I magnitude for UMi (left panel) and Dra (right panel). 
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In the lower panel of Fig. 6 the HB ridge line of UMi 
has been shifted to match the HB of Draco, that has been 
plotted with the same symbols used above. The fit is made 
difficult by the different HB morphology, but the solution 
obtained here with a shift of +0.42 mag in V and +0.02 
in V-I appears as the best possible match. The derived 
relative distance modulus is <^(to-m) = 0-42 ± 0.10. The 
result is in good agreement with the difference of moduli 
that can be obtained from Table 2 of Mateo (1998), i.e. 
<5(m-M) =0.47 ±0.18. 

4.2. Comparison with template clusters 

An indirect estimate of < Vrr > and of the distance 
modulus can be obtained by finding a match with the HB 
of well studied template clusters of similar metallicity, for 
which a robust direct estimate of < Vrr > is available 
(see, e.g. Montegriffo et al. 1998, and references therein). 

In Fig. 7 we report the results of the comparison with the 
photometry of M 68 {[Fe/H] = -2.09) by Walker (1994). 
The HB of UMi (upper panel) and Draco (lower panel) are 
plotted with the same symbols as Fig. 6, while the stars of 
M 68, shifted to match the HB of UMi and Dra, are repre- 
sented by open circles. The solid line is the mean level of 
the RR Lyrae in M 68 (< Vrr >= 15.64 + 0.01, according 
to Walker 1994) after the application of the shift needed 
to provide the match. 

A good match is found between the HB of M 68 and 
UMi applying a shift of 4.24 mag in V to the stars of M 68. 
Using the above quoted < Vrh > of M 68, and propagat- 
ing all the uncertainties < Vrr > {UMi) = 19.88 ± 0.10 
is inferred, in excellent agreement with the estimate by 
N88. We derive also the level of the Zero Age Hori- 
zontal Branch {Vzahb, a quantity that is best suited 
for comparison with theoretical models) by applying the 
same AV to the Vzahb value of M 68 reported by Fer- 
raro et al. (1999, hereafter F99), Vzahb = 15.75 + 0.05. 
We obtain VzAHsiUMi) = 19.99 + 0.11. Finally, tak- 
ing into account the difference in reddening and adopt- 
ing (m - M)o(M68) = 15.11 + 0.1 from F99 we obtain 
the following estimate for the distance modulus of UMi: 
(m — M)o = 19.38 + 0.16, inclusive of all the uncertainties. 

The HB morphology of Draco has no counterpart among 
the galactic globular clusters of similar metallicity. Hence, 
the fit with the HB of M 68 - shown in the lower panel of 
Fig. 7 -is worse than that found for UMi and do not provide 
a well constrained solution. However, with the adopted 
shift we obtain results that are in full agreement with the 
differential distance between UMi and Draco derived in 
§4.1. With the same procedure described above we obtain 
for Draco: < Vrr >= 20.30 + 0.12, Vzahb = 20.41 + 0.13 
and (m - M)o = 19.80 + 0.18 

In Fig. 8 are reported the results of the comparison with 
the HST photometry of M 92 {[Fe/H] = -2.24) by Ferraro 
and co-workers (Ferraro et al., in preparation), that is also 
in good agreement with the well calibrated ground-based 
photometry by Johnson & Bolte (1998). The symbols and 
the procedure are the same adopted in Fig. 7 for the com- 
parison with M 68. 

Also in this case a good match is obtained for UMi 
(upper panel), with the shifts reported in the figure. 
We adopt larger uncertainties in the derived quantities, 
with respect to the case of M 68, since M 92 is slightly 



more metal poor than UMi (and Draco). Adopting < 
Vrr > {M 92) = 15.08 + 0.01 from Kopacki (2001) we 
obtain < Vrr > {UMi) = 19.84 + 0.12, again in ex- 
cellent agreement with N88 as well as with the results 
from Fig. 7. Adopting Vzahb{M 92) = 15.30 + 10 
from F99 it is obtained VzAHB{UMi) = 20.06 + 0.14, 
in agreement with what found above. Finally, adopting 
(m - M)o(M 92) = 14.74 + 15 from F99, it is obtained 
(m - M)Q{UMi) = 19.47 + 0.20. 

M 92 appears to provide a better match to the HB 
of Draco, with respect to M 68 (Fig. 8, lower panel). 
With the same assuptions as above, we obtain for Draco: 
< Vrr >= 20.26 + 0.12, Vzahb = 20.48 + 0.14 and 
(m - M)o{Dra) = 19.89 + 0.20. The overall agreement 
with the results obtained from the comparison with M 68 
is fully satisfying. 

In conclusion, we adopt the weighted mean of the re- 
sults obtained above for < Vrr >, Vzahb and (m — M)o. 
The final results are < Vrr >== 19.86 + 0.09, Vzahb = 
20.02 + 0.09, for UMi, and < Vrr 20.28 + 0.10, 

Vzahb = 20.44 + 0.10, for Draco. The adopted dis- 
tance moduh are (m - M)o{UMi) = 19.41 + 0.12 and 
(m - M)o{Dra) = 19.84 + 0.14 

4.3. A consistency check: the RGB Tip 

The use of the Tip of the Red Giant Branch (TRGB) as 

a standard candle is a powerful technique to estimate the 
distances to galaxies that host old stellar populations (see 
Madore & Freedman 1998, and refereces therein). Bel- 
lazzini, Ferraro & Pancino (2001a) have recently provided 
a new robust calibration of the zero-point of the relation 
between the absolute I magnitude of the tip {Mj^'^^) 
and the metallicity that is independent of the distance 
scale based upon classical standard candles (RR Lyrae, 
Cepheids). This new calibration prompted us to start a 
large observational programme aimed at a re-assessment of 
the distances to Local Group galaxies based on the TRGB 
technique. The present study is the first of a series describ- 
ing the results of this program. Unfortunately, despite the 
large field sampled, UMi and Draco have too few stars (and 
a too low siirface brightness) to allow an individual safe 
application of the TRGB technique (see Sakai, Madore 
& Freedman 1996; Madore & Freedman 1995, hereafter 
MF95). However, it is of great interest for the final scope 
of the project to check if the distance of these galaxies is 
consistent with the new TRGB distance scale. 

To do this check we adopt the shift derived in §4.1 to 
move Draco at the distance of UMi and we try to detect 
the TRGB from the merged sample, having now a num- 
ber of stars nearly double with respect to each individual 
case. In the left panel of Fig. 9 we report the composite 
(I;V-I) CMD of the two galaxies together, after the appli- 
cation of the quoted shift to Draco. The right panels of 
Fig. 9 report the actual detection of the TRGB using the 
standard technique, as defined by Sakai, Madore & Freed- 
man (1996). The sharp cut-off of the Luminosity Func- 
tion (LF) of the RGB (that is the actual observational 
marker of the TRGB) is clearly visible in the histograms 
shown in the upper and middle right panels. The edge- 
detector Sobel filter shows the TRGB as the most pro- 
nounced peak at 1=15.57 + 0.05. The considered sample 
has 100 stars in the upper 1-magnitude bin, however, from 
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Fig. 7. — Matching the HBs of UMi (upper panel) and Draco (lower panel) with the HB of M 68. The symbols for UMi and Draco are the 
same as in Fig. 6 while the M 68 stars are plotted as open circles. The solid line is the average RR Ly level of M 68 sfithed by the amount 
specified in the upper left corners of the panels. 
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Fig. 8. — Matching the HBs of UMi (upper panel) and Draco (lower panel) with the HB of M 92. The symbols for UMi and Draco are the 
same as in Fig. 6 and 7 while the M 92 stars are plotted as open circles. The solid line is the average RR Ly level of M 92 sfithed by the 
amount specified in the upper left corners of the panels. 
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TRGB: UMi + Dra 




Fig. 9. — Left panel: composite CMD of UMi and Draco, obtained by shifting the magnitudes of Draco stars by -0.42 mag. Right panels: 
ordinary histogram (upper panel), generalized histogram (middle panel) and Sobel's filter response (lower panel) for the LF of the upper RGB 
of the composite sample. The position of the detected TRGB is indicated in the lower panel. 
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the CMD of field stars we estimate that ~ 20 % of them 
are likely foreground sources. Hence the actual number of 
stars in the uppermost 1-mag bin is in fact A''^ ~ 80. By 
means of numerical simulations MF95 showed that a rea- 
sonably safe detection of the TRGB can be made only if 
iV^ > 100, thus it is likely that our detection overestimates 
the actual I magnitude of the tip. To quantify this sys- 
tematic error we repeat the MF95 numerical experiment. 
We adopted the ridge line of M 68 as an RGB template 
and a model RGB-LF of the same form as the one used 
by MF95, imposing a cut-off at M/ = -4.00. Then we 
randomly generated 100 synthetic RGBs for each of the 
following cases: AT* = 20, 40, 60, 80, 100, 150, 200, and we 
detected the tip of each synthetic RGB with the standard 
technique. Observational errors extracted from a gaussian 
distribution with ay , cri similar to those of the observed 
stars were added to each synthetic star before the TRGB 
detection. The average Mi{TRGB) of the 100 simulated 
RGBs and the corresponding standard deviations are plot- 
ted as a function of A^* in Fig. 10. The results by MF95 are 
confirmed: for iV* > 100 the estimate from the standard 
technique is within < 0.02 mag from the true TRGB and 
the standard deviation is < 0.04 mag. The efficiency of the 
method drops suddenly for iV^, < 80 and the associated un- 
certainty becomes ~ 0.1 or larger. The difference between 
the true luminosity of the tip and the observed one can be 
modeled with a first order polinomial as a function of N^,, 
for < 100: 

SItrgb = Itrgb - Itrgb = 0.00227V, - 0.236 

In the present case, with N-k ~ 80, SItrgb = —0.07 ± 
0.08, thus our final estimate for the TRGB location is 
Itrgb = 15.50 ± 0.14, where we have considered also the 
uncertainties associated with the correction and with the 
Dra-UMi shift we have applied to obtain the considered 
merged sample. The galaxies are dominated by a metal 
poor population quite similar to that of the template used 
by Bellazzini, Ferraro & Pancino (2001a), that obtained 
MJRGb ^ _4 04 ± 0.12 at [Fe/H] = -1.7, for the globu- 
lar cluster w Centauri. Adopting E{B ~V) = 0.03, we ob- 
tain (m-M)o = 19.50±0.20 and (m-M)o = 19.92±0.27 
for UMi and Draco respectively, in good agreement with 
the results of §4.2. 

Hence we conclude that the distance moduli derived in 
§4.2 are fully consistent with the TRGB distance scale as 
calibrated by Bellazzini, Ferraro & Pancino (2001a). 

5. THE STELLAR CONTENT OF UMI AND DRACO 

5.1. The RGB hump 

The RGB-Bump is an evolutionary feature occurring 
along the RGB which flags the point where the H-burning 
shell crosses the chemical discontinuity left by the max- 
imum penetration of the convective envelope. Theoreti- 
cal models predict that the position of the RGB-Bump is 
mainly driven by metallicity, with a mild dependence on 
age (see F99 and references therein). 

From the observational point of view the feature was 
observed for the first time in a globular cluster by King, 
Da Costa & Demarque (1985) and it was only recently 
identified in a significant number of clusters (Fusi Pecci et 
al. 1990; Zoccali et al. 1999, F99). The first detection of 



the RGB-Bump in a galaxy is very recent too. Majewski 
et al. (1999) identified a double RGB-bump in the Sculp- 
tor dwarf spheroidal and took this evidence as indicative 
of the presence of two populations of different metallici- 
ties. A very similar feature has been revealed also in the 
Sextans dSph by Bellazzini, Ferraro & Pancino (2001b). 
At present, these are the only two detections of the RGB- 
bump ever obtained in stellar systems other than globu- 
lars. 

The change in the slope of the cumulative luminosity 
function (LF) and the excess of star counts in the differ- 
ential LF of the RGB are the main tools to identify the 
RGB-Bump. In particular Fusi Pecci et al. (1990) suggest 
that the change in the slope of the cumulative LF is the 
safest indicator of the RGB-Bmnp. In Fig. 11 the differ- 
ential and cumulative LFs of the RGB in the region of the 
bump are shown for both galaxies. 

The RGB-bump is clearly detected in UMi at Vf'"™^ = 
19.40 ± 0.06. This is the first single RGB-bump ever de- 
tected in a galaxy. On the other hand, there is no sign of 
this feature in the LF of Draco. We argue that a larger 
spread of metallity and/or age in the stellar population of 
Draco, with respect to UMi, is responsible for the smearing 
of the RGB-bump feature in this galaxy (see §5.2). This 
suggests that these dSphs - that are quite similar under 
many aspects (see Bellazzini, Ferraro & Pancino 2001b, 
and references therein) - had different evolutionary histo- 
ries at early times. 

The magnitude difference between the RGB-bump and 
the HB level (AU^™^) has been used by Fusi Pecci et 
al. (1990), and later by F99 and Zoccah et al. (1999), to 
compare the observations with the theoretical predictions. 
The AV^'^^ parameter is mainly a function of the metal- 
licity and it has only a mild dependence on age. Adopt- 
ing the ZAHB levels obtained in §4.2, we obtain for UMi 
AV^j^^ = —0.62 ± 0.11. Since it has been shown that 
UMi is dominated by a very old population that formed 
in a rather short timescale (Carrera et al. 2002; Mighell 
& Burke 1999; Feltzing, Gilmore & Wyse 1999; Dolphin 
2002) we can safely use the AV^j^^ parameter as an indi- 
cator of the mean metal content of the old stars. We use 
the data by F99 to obtain a fit of [Fe /H] as a fimction of 
AV^j^^, adopting the metallicity scales by Zinn & West 
(1984) (ZW) and Carretta & Gratton (1997) (CG). It is 
found: 

[Fe/H]zw = 1.426AV^"™P - 1.233 {rms = 0.08) 

[Fe/H]cG = -0.430AV^7*f%1.183AV^"^f-1.046 {rms = 0.07 

With these relations we obtain for UMi < 
[Fe/H]zw >= "2.1 ± 0.2 and < [Fe/H]cG >= 
— 1.95 ± 0.2, in good agreement with the estimates ob- 
tained with other methods (see Mateo 1998; Carrera et al. 
2002; Shetrone, Cote & Sargent 2001). 

5.2. The metallieity distributions 

In Fig. 12 the (I,V-I) distribution of the upper RGB 
of UMi (left panel) and Draco (right panel) are com- 
pared with the ridge lines of the template globular clusters 
NGC 6341, NGC 6205 and NGC 288, from left to right. 
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Fig. 10. — Detected absolute I magnitude of the TRGB as a function of the number of stars in the upper 1 mag bin. Each point is the 
average of 100 simulated samples having Af* stars in the upper bin. The error bars are the corresponding standard deviations. 
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Fig. 11. — DifTcrcntial (upper panels) and cumulative (lower panels) LF of the RGB of UMi (left panels) and Draco (right panels). The 
position of the RGB-bump of UMi is clearly indicated. 
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taken from the homogeneous set by Saviane et ah (2000). 
In the Zinn & West (1984) metaUicity scale the templates 
have [Fe/H] = -2.24, -1.65 and -1.40 respectively, while 
in the Carretta & Gratton (1997) scale their metallicities 
are [Fe/H] = -2.16,-1.39 and -1.07 respectively. The 
metallicities, distance moduli and reddenings of the tem- 
plates are taken from F99. For the two galaxies we adopted 
the estimates of §4. The meaning of the horizontal lines 
and of the large empty symbol will be discussed in §5.2.1. 

The large majority of the UMi RGB stars are enclosed 
between the ridge lines of NGC 6341 and NGC 6205, im- 
plying a mean metaUicity around [FelH]zw — —2.0 in 
good agreement with the result found in §5.1 as well as 
with previous estimates, either photometric or spectro- 
scopic (sec Carrcra ct al. 2002; Shetrone, Cote & Sargent 
2001, and references therein). The RGB of Draco seems 
slightly redder than the one of UMi, suggesting a mean 
metaUicity [Fe/ H]zw — ^l-^, also in agreement with 
previous estimates (see Aparicio, Carrcra & Martinez- 
Delgado 2001; Shetrone, Cote & Sargent 2001, and ref- 
erences therein). Both galaxies result significantly more 
metal rich (by ~ 0.3 dex) if the CG metaUicity scale is 
adopted, as a natural consequence of the higher metaUicity 
value of all the adopted templates in this scale. Indepen- 
dently of the assumed scale, it is clear that the color spread 
along the RGB of both galaxies is significantly larger than 
the photometric scatter (see Fig. 4), thus we confirm the 
presence of an intrinsic metaUicity spread (Shetrone, Cote 
& Sargent 2001; Carrera et al. 2002; Aparicio, Carrera & 
Martinez-Delgado 2001). 

The RGB color of old populations is mainly driven by 
the metal content of the stars. Based on this principle, 
photometric Aletallicity Distributions (MD) can be ob- 
tained by suitable interpolation between a grid of metal- 
licity templates, a technique that has been widely applied 
in recent years (see, e.g. Holland, Fahlman & Richer 1996; 
Harris & Harris 2000; Saviane et al. 2000, and references 
therein). Here we obtained the MDs shown in Fig. 13 
and Fig. 14 by interpolating in color between the RGB 
ridge-lines of the adopted templates (see Fig. 12) for the 
stars in the magnitude range —2.9 < Mj < —3.9. While 
the accuracy of the single metaUicity estimate is low, the 
main characteristics of the distribution as a whole (as, for 
instance, its mean and dispersion) are quite robust, being 
based on a large number of stars. This technique is partic- 
ularly appropriate in the present case since we know that 
both galaxies are dominated by very old stellar popula- 
tions (Mighell & Burke 1999; Grillmair et al. 1998; Dolphin 
2002; Aparicio, Carrera & Martinez-Delgado 2001; Car- 
rera et al. 2002). However, the actual metaUicity spread 
dcaived with this tc;c;lmique (a) is the convolution of the 
intrinsic metaUicity spread with the color spread due to ob- 
servational scatter and, (b) the impact of this latter factor 
depends on the mean metaUicity, since the same observa- 
tional scatter should produce a larger (apparent) metaUic- 
ity spread in the low metaUicity regime, where the RGB 
ridge lines are steeper, with respect to the high-metallicity 
regime. In fact, A[Fe/H]/ A{V—I) decreases with increas- 
ing metaUicity. To deal with this problems we compare the 
observed photometric MDs with those of the recovered ar- 
tificial stars shown in Fig. 4, that closely model the CMD 
of a SSP observed and reduced under the same conditions 



as real stars. 

In Fig. 13 the MD of Draco (right panels) and the MD 
of corresponding synthetic SSP (left panels; see Fig. 5) are 
compared. All the MDs are normalized to their maximum 
value. The upper panels show the histograms in the ZW 
scale. The total number of included stars is reported in the 
upper right corner. The middle and lower panels show the 
generalized histograms (continuous lines) in the ZW scale 
and in the CG scale, respectively (see Holland, Fahlman 
& Richer 1996; Laird et al. 1988, and references therein). 
The measured standard deviation (cr) is reported in the 
upper right corner of each panel. The MD of the synthetic 
SSP represents the response of the whole procedure to the 
pure observational scatter, or in other words is the ob- 
served version of a true metaUicity distribution having the 
form of a Dirac delta function. The plots in the left panels 
of Fig. 13 demonstrate that the observational scatter pro- 
duce a sizeable (spurious) metaUicity spread. To obtain 
a sensible measure of the true intrinsic metaUicity spread 
we have to deconvolve this "response function" from the 
observed distributions shown in the right panels of Fig. 13. 
The dotted lines are the gaussian distributions obtained by 
deconvolving the response function from a gaussian dis- 
tribution having the same mean and standard deviation 
of the observed distributions. The standard deviation of 
the dotted curves is a good estimate of the true intrinsic 
metaUicity spread ai. 

In the ZW scale the MD of Draco spans more than 0.6 
dex. The mean and median metallicities are < [Fe/H] >= 
-1.7±0.1, while the mode is [Fe/H]mod = -1.6±0.1. The 
intrinsic 1 — cr scatter is CTj = 0.13 dex. There is a marginal 
indication of a secondary peak at lower metallicty, but 
the statistical significance of this feature is low, thus we 
don't comment further on it. These estimates are in good 
agreement with the results of the spectroscopic analysis 
by Lehnert et al. (1992) and Shetrone, Cote & Sargent 
(2001). The detailed comparison between our photomet- 
ric estimates and the study of Shetrone, Cote & Sargent 
(2001) is discussed in §5.2.1, below. Our method is insensi- 
tive to very metal poor stars {[Fe/H] < —2.5), so they are 
obviously not represented in our metaUicity distribution. 
However from Fig. 12 it may be evinced that the number 
of possible RGB stars bluer (i.e. more metal poor) than 
the NGC 6341 template (at [Fe/H] = -2.24) should be 
small in both galaxy (see §5.2.1 for further details). 

In Fig. 14 the MDs of UMi and of its synthetic SSP are 
reported in the same way as in Fig. 13. It can be appre- 
ciated that though the photometric scatter is very similar 
for Draco and UMi (see Fig. 4) the response function has 
a larger a in the case of UMi. This is due to larger sensi- 
tivity of metaUicity to color at the most metal poor regime 
described above. In the ZW scale the MD of UMi spans 
a range nearly as wide as that of Draco. The mean and 
median metallicities are < [Fe/H] >= —1.8 ± 0.1, while 
the mode is [Fe/H]mod — — 1.9±0.1. The intrinsic scatter 
is CTj = 0.10 dex. Also in this case the agreement with 
results of other studies (Shetrone, Cote & Sargent 2001; 
Aparicio, Carrera & Martinez-Delgado 2001) and with the 
constraints obtained in §5.1 is quite good. 

From the comparison of the MDs of the two galaxies we 
can conclude that Draco is slightly less metal deficient than 
UMi, in average, and it presents a slightly larger metallic- 
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Fig. 12. — The (I,V-I) distribution of the upper RGB of UMi (left panel) and Draco (right panel) are compared with the ridge lines of 
the templates globular clusters NGC 6341 {[Fe/H]zw = -2.24; [Fe/H]cG = -2.16), NGC 6205 {[Fe/H]zw = -1-65; [Fe/H]cG = -1.39) 
and NGC 288 {{Fe/H]zw = -1-40; [Fe/H]cG = -1-07), from left to right. The long-dashed line is the expected level of the TRGB with 
the adopted distance moduli (from §4.2) and the TRGB calibration by Bellaazini, Ferraro & Pancino (2001a). The dotted lines enclose the 
full range of uncertainty in the expected magnitude of the TRGB. The large empty symbols marks the stars observed by Shetrone, Cote 
& Sargent (2001) acoording to the metaJlicity estimated by these authors. Circles: {Fe/H\ < —2.10, pentagons: —2.10 < [Fe/H] < —1.8, 
squares: -1.8 < [Fe/H] < -1.5, triangles: -1.5 < [Fe/H] < -1.3. 
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Fig. 13. — The photometric metallicity distribution of Dreico (right panels) are compared with the distribution obtained with the same 
technique from the synthetic SSP we adopted for the artificial stars experiments (left panels). Upper panels: odinary histograms in the 
ZW metallicity scale. Middle panels: generalized histograms in the ZW metallicity scale. Lower panels: generalized histograms in the CG 
metallicity scale. The adopted smoothing lenght of the generalized histograms is equal to the measured a of the response function, i.e. the 
natural smoothing lenght. 
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Fig. 14. — The photometric metallicity distribution of UMi (right panels) are compared with the distribution obtained with the same 
technique from the synthetic SSP we adopted for the artificial stars experiments. The symbols are the same as in Fig. 13. 
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ity spread. The parameters obtained from the photometric 
MDs (for Draco and UMi and for both metalhcity scales) 
are summarized in Table 3. 

The MDs in the CG scale are remarkably similar to the 
MDs in the ZW scale described above, but, as expected, 
arc shifted to higher metallicitics and show a sUghtly larger 
spread, thus the conclusions drawn above are unchanged. 
On the other hand, it is important to remark that the 
comparison of the observed distributions with chemical 
evolution models would lead to different conclusions de- 
pending on the assumed metalhcity scale. In particular, 
by adopting the CG scale, one should conclude that the 
first stars in Draco and UMi formed from an ISM already 
relatively enriched and reached a more advanced stage of 
chemical evolution than suggested by the ZW scale. This 
kind of uncertainties in the metalhcity scales may have 
other dangerous drawbacks, hampering our interpretation 
of astrophysical data in several cases (see, e.g. Catelan et 
al. 2001). 

5.2.1. Comparison with spectroscopic estimates 

A detailed comparison with the recent results by 
Shetrone, Cote & Sargent (2001) may be useful to 
check the consistency between high-resolution spectro- 
scopic abundances and photometric metalhcity estimates. 
The stars observed by Shetrone, Cote & Sargent (2001) 
arc plotted in Fig. 12 as large empty symbols. Different 
symbols arc adopted according to the metalhcity estimates 
by Shetrone, Cote & Sargent (2001). Circles are stars 
with [Fe/H] < -2.10, pentagons have -2.10 < [Fe/H] < 

— 1.8, squares have —1.8 < [Fe/H] < —1.5, and trian- 
gles have —1.5 < [Fe/H] < —1.3 (actually, there are only 
two stars in this latter bin both having [Fe/H] — —1.45). 
In both panels of Fig. 12, the dashed line is the expected 
level of the TRGB with the adopted distance moduli (from 
§4.2) and the TRGB calibration by Bellazzini, Ferraro & 
Pancino (2001a). The dotted lines enclose the full range 
of uncertainty in the expected magnitude of the TRGB, 
including observational and calibration uncertainties. 

First of all we note that all the considered stars have 
I ^ Itrgb, to within the uncertainties. This is com- 
patible with the hypothesis that they are first ascent 
RGB stars, as suggested by Shetrone, Cote & Sargent 
(2001). Second, the metalhcity rank derived from spec- 
troscopy is correctly reproduced by our photometry: stars 
with higher spectroscopic metalhcity estimates lie on red- 
der RGB loci. Finally, the zero-point of the spectro- 
scopic scale is in good agreement (within the uncertain- 
ties) with that derived from our photometry, if the CG 
metalhcity scale is adopted. Three out of four of the 
stars in the more metal poor metalhcity bin (circles) have 
-2.36 ± 0.09 < [Fe/H] < -2.17 ± 0.12 and lie around 
the ridge line of NGC 6341 {[Fe/H]cG = -2.16 ± 0.10). 
The two more metal rich stars (triangles, [Fe/H] = 

— 1.45 ± 0.07) are slightly bluer than the ridge line of 
NGC 6205 {[Fe/H]cG = -1.39 ± 0.10). Hence both the 
metal rich and the metal poor zero-points seems to agree 
quite well. The only exception is the star ^ 119 in Draco 
that is slightly redder than the ridge line of NGC 6341 
and that, according to Shetrone, Cote & Sargent (2001) 
has [Fe/H] = — 2.97±0.15. We have no suggestion to 
reconcile such a large discrepancy between color and spec- 



troscopic metalhcity, except for trying to repeat both mea- 
sures. 

If the ZW scale is adopted for the templates shown in 
Fig. 12 the consistency in the metal poor zero-point is 
preserved {[Fe/H]zw = -2.24 ± 0.20 for NGC 6341) but 
a serious problem emerges at the metal rich end, where 
the two stars with [Fe/H] = —1.45 lie to the blue of the 
[Fe/H] = -1.65 template and the [Fe/H] = -1.4 tem- 
plate (NGC 288) is more than 0.1 mag redder than these 
stars. The problem with star =ff 119 is not mitigated by 
the assumption of the ZW metalhcity scale. Thus it ap- 
pear that the adoption of the GC metalhcity scale pro- 
vide a much higher degree of consistensy between spec- 
troscopic and photometric metalhcity estimates, at least 
in the present case. The above comparison may suggest 
that the metallicitics by Shetrone, Cote & Sargent (2001) 
are consistent with the CG scale which is also based on 
high-resolution measures. 

Finally, another possible source of uncertainty in the 
photometric metalhcity estimates is the small difference 
in the average abundance of a-elements between Dra and 
UMi and the template clusters (see Shetrone, Cote & Sar- 
gent 2001, and discussion therein). 

5.3. Populations gradients 

Radial population gradients are ubiquitous in dwarf 
galaxies (see Harbeck et al. 2001; Saviane et al. 2001; Hopp 
et al. 2000; Tosi et al. 2001, and references therein). Al- 
most in all cases old/metal poor stars are found to be more 
abundant in the outer regions of dwarf galaxies, while more 
metal rich and/or younger stars are preferentially found 
near center of galaxies (Harbeck et al. 2001). Here we 
look for population gradients by using the fraction of blue 
HB stars (i.e. those with {V — I) < 0.4) in different con- 
centric radial annuli. The small number of available HB 
stars limits the spatial resolution to just three annuli. The 
results are shown in Fig. 15. In UMi (upper panel), the 
fraction of BHB stars appears to rise slightly in the in- 
ner 4 arcmin. However the effect is at less than la, and 
the data are consistent with no population gradient, in 
agreement with what found by Carrera et al. (2002). The 
opposite effect is found in Draco (lower panel): the BHB 
fraction drops by a factor 2 in the inner 4 arcmin, i.e. the 
classical gradient described above, with more metal rich 
and/or younger stars clustered in the central region. The 
significance of the result is just marginal (2a-), thus not in 
disagreement with the null detection reported by Aparicio, 
Carrera & Martinez-Delgado (2001). 

6. STRUCTURE 

We take advantage of our well defined CMDs to select 
likely members of UMi and Draco with which reliable den- 
sity contour maps can be obtained. The upper panels of 
Fig. 16 show the adopted selections, while in the lower pan- 
els of the same figure the isodensity contour maps of the 
selected members of UMi and Draco are shown. Our maps 
cover the inner regions of the galaxies, within 1-2 core radii 
(jc), and are obtained by computing the stellar density on 
an uniform grid with points spaced by 100 px. At each 
point of the grid the density is computed over a circle of 
radius r = 3 arcmin. This choice provides a suitable degree 
of smoothing and minimizes the effects of possible small 
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empty regions (f.i, saturated stars or dead CCD columns). 
The continuous hues are isodensity contours starting from 
1 stars /arcmin^ and spaced by 0.5 stars /arcmin^. The 
innermost contour corresponds to 4.5 stars /arcmin^ for 
UMi and to 6.0 stars/ arcmin^ for Draco. To give a term 
of comparison Klcyna ct al. (1998) - using a sample with 
similar limiting magnitude - found a very similar value 
for the maximum density of UMi (~ 4 stars /arcmiri^) 
and estimated the density of the background to be 0.31 
stars I arcmir? at a distance of ^ 30 arcmin from the cen- 
ter of UMi, much farther out than our outermost isodense. 

The structure of the two galaxies, as shown in Fig. 16, 
is in excellent agreement with previous results (Irwin & 
Hatzidimitriou 1995; Kleyna et al. 1998; Piatek et al. 2002; 
Odenkirchcn ct al. 2001): the cllipticity, position angle and 
the ratio between the central surface brightness are the 
same as reported in the review by Mateo (1998). Draco is 
more dense and has a much smoother profile with respect 
to UMi. Its isodensity contours are quite regular, show- 
ing just a marginal trend toward more boxy shapes in the 
outer parts. On the other hand the inner contours of UMi 
appear quite structured and asymmetrical and deserve a 
deeper analysis and discussion. 

The detection of small scale clustering and asymmet- 
ric structures in the inner parts of UMi has been claimed 
by many different authors since the early '80 (see, e.g. 
Olszewski & Aronson 1985; Irwin & Hatzidimitriou 1995; 
Klcyna ct al. 1998; BattincUi & Demers 1999; Dcmcrs & 
Battinelh 2001; Eskridge & Schweitzer 2001), but a clear- 
cut demonstration of the statistical significance of such 
structures is still lacking. In Fig. 16 the isodensity con- 
tours up to 3.0 stars jar emir? appear quite elongated but 
similar in shape and symmetrical. On the other hand the 
peak of density {X ~ 1.8;y ~ 3) is clearly off-centered 
with respect to the center of symmetry of the outer con- 
tours (X ~ 4.3;y ~ 5.3). We note that the density con- 
trast between the density peak and the last symmetric 
contour is larger than 4 cr, i.e. quite significant. It is also 
remarkable that the off-centered main density peak dis- 
plays a much rounder shape, with respect to the rest of 
the galaxy. On the other hand, the dotted contour (corre- 
sponding to a density of 3.2 star s / ar emir?) is reported to 
show that the substructures in the NE region of the galaxy 
that are clearly visible in the maps by Irwin & Hatzidim- 
itriou (1995) and Kleyna et al. (1998) are present also in 
our sample but have a low density contrast with respect 
to the last symmetric contour, thus low statistical signifi- 
cance, as concluded also by Kleyna et al. (1998). 

It is important to remark that the only other spheroidal 
galaxy of the Local Group that has a clearly off-centered 
density peak is the Sagittarius dSph, whose structure is 
greatly strained by the Galactic tidal field (see Ibata et 
al. 1997; Newberg et al. 2002). Independent evidence of 
the effects of the Galactic tidal fields on UMi has been 
reported by Martinez-Dclgado et al. (2001). 

Kleyna et al. (1998) found significant asymmetry in the 
distribution of stars along the major axis. Here we con- 
centrate on the clustering properties of UMi stars. We de- 
fine dn as the distance of a given star to its n—th nearest 
neighbour. In the following we will show the results ob- 
tained with ^200) i-e. the distance to the 200— th nearest 
neighbour, but we remark that the results are preserved 



for a very wide range of n (50 < n < 500). In the up- 
per panel of Fig. 17 the distributions of d2oo are shown 
for three template samples having the same dimensions 
of the observed ones (~ 2500 member stars per galaxy). 
The dotted curve is the (i2oo distribution of a sample ran- 
domly drawn from a uniform distribution. The continuous 
line and dashed lines are the c?2oo distributions of sam- 
ples drawn from an elliptical Gaussian distribution having 
= rdUMi) and (Jx/<7y = {a/b)uMi, and = rdDra) 
and (Tx/cTy = {a/b)]jra, respectively, where a and b are 
the semi-major and semi-minor axis of the galaxies, as 
listed by Mateo (1998). It can be easily appreciated that 
(1) more concentrated (denser) systems have ^200 distri- 
butions peaked at shorter lenghts and (2) independently 
of the assumed density profile, the fi2oo distributions of 
symmetric and smooth systems show a single major peak, 
marking the characteristic clustering scale, followed by a 
very sharp cut-off. In these systems there are no stars 
having d2oo shorter than the cut-off value. 

In the middle panel of Fig. 17 the d2oo distribution of 
Draco is presented. The observed distribution is very well 
reproduced by its template: the characteristic clustering 
scale is ~ 500 px ~ 2.3' and no star with clustering scale 
shorter than the cut-off (~ 480 px) is observed. On the 
other hand, the d2oo distribution of UMi (shown in the 
lower panel of Fig. 17) is quite different from that of Draco 
and from the smooth/symmetric templates. The major 
peak and the cut-off arc well reproduced by the corre- 
sponding template but there is a clear secondary peak at 
shorter scale lenghts with respect to the main cut-off. The 
presence of this secondary peak shows that there are two 
characteristic scales of clustering in UMi, a longer one 
(2± 700 px ~ 3.2') associated with the general properties of 
the system and a shorter one (~ 600 px ~ 2.7') associated 
with the more compact structures corresponding, to the 
off-centered density peak. 

Can the observed d2QQ distribution of UMi arise by 
chance (i.e. small number fluctuations) from an intrin- 
sically smooth and unstructured system ? To answer 
this question we extracted 1000 random samples from the 
model assumed for the UMi template shown in the upper 
panel of Fig. 17 and for each of them we computed the 
fraction of stars having ^200 lower than the observed cut- 
off {F{d2oo < ^'200"°^''^' hereafter Fd, for brevity), fixed at 
'^200 = 630 px. The distribution of the Fd of the simu- 
lated samples is shown in Fig. 18. The observed fraction 
{Fd = 0.12) is indicated by a dotted line. Fig. 18 shows 
that (a) in no case a sample having Fd equal or greater 
than the observed value has been extracted from the un- 
structured model, (b) the highest Fd value reached by a 
simulation {Fd ~ 0.08) is significantly lower than the ob- 
served one, and (c) 95 % of the simulated samples have 
Fd < 0.04. We conclude that the observed presence of two 
characteristic clustering scales is a real and statistically 
significant property of UMi. 

Finally, we find that the spread in magnitude of the blue 
HBs of Draco and UMi are small and compatible with the 
observational scatter at that magnitude. Thus we can't 
sec any sign of the significant elongation along the line of 
sight that was required in the model by Klessen & Kroupa 
(1998) to explane the large velocity dispersion observed 
in dSph galaxies without the need of massive dark mat- 
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Fig. 16. — Upper panels: CMDs of the two galaxies. The stars selected as members (bold face dots) are used to derive the density contour 
maps shown in the lower panels. The outermost contour corresponds to a density of 1 star arcmin"^; the contours are plotted in steps of 
0.5 stars arcmin"^ (i.e. 1.0, 1.5, 2.0, etc.). The innermost contour corresponds to a density of 4.5 stars arcmin"^ for UMi, and 6.0 stars 
arcmin"^ for Draco. The dotted contour in the map of UMi corresponds to a density of 3.2 stars arcmin"^. It has been reported to show 
the low contrast structures in the north-eastern region of the galaxy that have already been noted by Irwin & Hatzidimitriou (1995). The 
dotted contour shows also that the symmetry seen in the outer isodense curves brealjs just above the density of 3 arcmin"^. The density is 
estimated over a 100 pixel spaced grid. At each point of the grid the density is estimated in a circle with radius = 3 arcmin. 
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ter halos (see also Aparicio, Carrera & Martinez-Delgado 
2001). 

7. DISCUSSION 

The distance estimates to Draco and UMi derived in §4 
are tied to the distance scale introduced by F99 and com- 
patible with the TRGB scale introduced by Bellazzini, Fer- 
raro & Pancino (2G01a). The F99 scale is fully consistent 
with the scale based on the revised Hipparcos parallaxes 
by Carretta et al. (2000). F99 showed that their distance 
moduli of globular clusters arc tipically ~ 0.2 mag larger 
than those tied to traditional Mv{RR) — [Fe/H] calibrat- 
ing relation. This is probably the main reason why we 
find distance moduli for UMi and Draco that are larger 
by 0.2 — 0.3 mag than what generally found in literature 
(see, e.g. Mateo 1998, and references therein). We note 
also that, independently of the distance scale, our proce- 
dure of HB matching leads to results that are in excellent 
agreement with works based on the RR Lyrae (e.g., N88 
Nemec 1985) while a sensible mismatch is noticed with re- 
spect to the Main Sequence fitting technique adopted by 
MigheU & Burke (1999), using HST data. We argue that 
the problems with the absolute photometric calibration of 
HST-WFPC2 data may be at the origin of the observed 
differences. 

The comparison with the recent distance estimates by 
Aparicio, Carrera & Martinez-Delgado (2001) (for Draco) 
and Carrera et al. (2002) (for UMi) is of particular in- 
terest. Aparicio, Carrera & Martinez-Delgado (2001) de- 
rived a distance modulus to Draco by (1) estimating Vrr 
from the observed HB and (2) by adopting the classical 
Mv{RR) - [Fe/H] relation by Lee et al. (1999). They 
obtain Vrr = 20.14 ± 0.12, compatible with our estimate 
within the uncertainties, but their final distance modulus 
((m — M)o = 19.5 ± 0.2) is 0.34 mag smaller than ours. 
However, if the above quoted systematic difference of ~ 0.2 
mag between the F99 distance scale and the scale by Lee 
et al. (1999) is taken into account the consistency between 
the two estimates is clearly recovered. On the other hand, 
Carrera et al. (2002) derived the mean level of RR Lyrae 
in UMi by a comparison with template globular cluster of 
similar metallicity, i.e. the same method adopted here. 
Moreover, the distance moduli of the template clusters 
were taken from Reid (1999) who derived them from Main 
Sequence fitting to Hipparcos subdwarfs. With this ap- 
proach Carrera et al. (2002) obtained VrR = 19.84 ± 0.07 
and (m — M)o = 19.40 ± 0.10 in excellent agreement with 
our results. 

Independently of the above considerations, it is impor- 
tant to remark that the actual uncertainty on state-of-the- 
art distance estimates to Draco and UMi remains ~ ±25 
%, in spite of all efforts. It is clear that classical distance 
indicators are unefhcient in this context and valid alterna- 
tives are certainly needed. The search for detached double- 
lined eclipsing binaries (DD-EB) is now possible with new 
generation instruments and may be very rewarding, given 
the high binary fraction that has been reported for these 
two galaxies (Olszewski, Pryor & Armandroff 1996). 

7.1. Evolutionary history 

The metallicity distributions presented here confirm 
(with larger samples, with respect to existing spectroscopic 



surveys) that the stellar populations of both galaxies have 
a sizeable spread in metal content. Since they are dom- 
inated by very old stars, we shall conclude that signifi- 
cant self-enrichment took place in these systems in a quite 
short time scale (i.e., few Gyr) at very early epochs. The 
same conclusions have been previously reached for Sculp- 
tor (Majewski et al. 1999) and Sextans (Bellazzini, Ferraro 
& Pancino 2001b). Shetrone, Cote & Sargent (2001) have 
shown that some of the stars of UMi, Dra and Sex show 
abundance patterns significantly different from those typ- 
ical of Galactic halo stars. In particular it has been sug- 
gested (Ikuta & Arimoto 2002) that the transition from an 
a -enhanced regime to nearly solar [a/Fe] ratios occurred 
at a much lower [Fe/H] value than in our own Galaxy. This 
evidence led Shetrone, Cote & Sargent (2001) to conclude 
that it is unlikely that systems like UMi and Draco gave 
a significant contribution to the assembly of the Galactic 
halo. We remark here that this conclusion is correct only 
for satellites that would have been accreted by the Milky 
Way in an advanced stage of their chemical evolution. A 
very early accretion (i.e., within ~ 1 Gyr from the on- 
set of star formation) or simply a very early stripping of 
the gaseous component by the Galaxy wouldn't leave any 
peculiar chemical signature in a Galactic halo largely com- 
posed by building blocks like Draco, UMi or Sextans (as 
they were in their first Gyr of "stellar" life). 

Our results show also that, while the average metallicity 
of the two galaxy is similar, the difference in the MDs is 
significant. The MD of UMi peaks at [Fe/H] c:i -1.9 and 
barely reaches [Fe/H] ~ — 1.6 while that of Draco peaks 
at [Fe/H] ~ -1.6 and reaches [Fe/H] ~ -1.3. This fact 
has to be taken into account, in particular in the interpre- 
tation of the great difference in HB morphology between 
the two galaxies. Dra and UMi are not a good "second 
parameter pair" (see Bellazzini et al. 2001, and references 
therein): a sizeable part of the HB morphology difference 
may be due to the differences in the actual MD. The large 
population of binaries (and Blue Stragglers) hosted by the 
galaxies (Olszewski, Pryor & Armandroff 1996) may also 
be have some impact on the observed HB morphology, es- 
pecially in the case of Draco (sec Fusi Pecci et al. 1992; 
Bellazzini et al. 2001; Bellazzini et al. 2002a; Carney et 
al. 2001, and references therein). 

7.2. A striking structural difference: not so twins after 

all? 

Draco and UMi have similar distance from the center 
of the Galaxy, similar luminous mass and similar star for- 
mation histories. Both are thought to have high M/L ra- 
tios and probably have also similar orbits, since they have 
similar radial velocities and (possibly) proper motions (see 
Mateo 1998; Scholz & Irwin 1994; Schweitzer, Cudworth 
& Majewski 1997). In spite of that, while Draco shows a 
very symmetrical and smooth profile out to a very large 
distance from the center (Odenkirchen et al. 2001; Piatek 
et al. 2002), we have demonstrated that even the innermost 
regions of UMi are strongly structured and asymmetrical. 
The presence of a massive and extended CDM halo should 
have inhibited the formation of any significant substriic- 
ture in the stellar component of UMi. The star formation 
activity cannot be responsible of the disturbed structure of 
this galaxy since significant star formation episodes have 
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ceased many Gyr ago. Furthermore we accurately checked 
for possible differences in the stellar populations between 
the off-centered (and nearly round) density peak of UMi 
and its outer regions and we found none. 

Ferraro, Bellazzini & Pancino (2002) showed that a 
small self-gravitating stellar system is embedded into the 
main body of the peculiar globular cluster w Cen. This 
stellar system has likely managed to preserve its individ- 
uality over many Gyr, while orbiting into the stellar halo 
of the cluster. A similar possibility could be envisaged 
also to explain the origin of the peculiar density peak of 
UMi. The process of hyerarchical merging is expected to 
be scale-free, thus the records of the formation of a galaxy 
like UMi from smaller fragments may be still detectable 
today (sec Fellhauor & Kroupa 2000), possibly favoured 
by the fact that the putative fragment we are considering 
is much denser than the stellar medium it is embedded in. 
The main argument against this hypothesis is the already 
reported observational evidence that the stellar popula- 
tion in the density peak docs not differ from the average 
UMi population (as opposite to the to Cen case, see Fer- 
raro, Bellazzini fc Pancino 2002, and references therein). 
Large high-resolution spectroscopic surveys may possibly 
provide more stringent indications on the viability of this 
scenario. 

While we are writing, a new preprint was posted (Palma 
et al. 2002) in which our results about the structure of UMi 
are fully confirmed'''. Furthermore Palma ct al. (2002) were 
able to follow the structure of UMi out to a distance of 
^ 1.7 deg from the center and demonstrated that (1) the 
galaxy is elongated along the direction of its proper motion 
vector and (2) the outer isodensity contours of UMi have 
the typical S-shape of tidally disturbed stellar systems. If 
it will be spectroscopically confirmed that the stars in the 
S-shaped structures are unbound, or the signature of ap- 
parent rotation will be found, this would directly disprove 



a specific prediction of the CDM simulations (Hayashi et 
al. 2002). In any case, the strongly disturbed structure of 
UMi provides a severe challenge to standard CDM scenar- 
ios, above all if the comparison with Draco is considered^. 
Moreover, Palma et al. (2002) showed that previous es- 
timates of the total luminosity of UMi had missed a sig- 
nificant amount of light that is found at large distances 
from the center of the galaxy. Their estimate of total lu- 
minosity is 2.7 times larger than previous ones. Taking 
this into account and considering the (possible) effects of 
anisotropic velocity dispersions they were able to reduce 
the mass to light ratio of UMi down to M/L ~ 16. If we 
include in this computation the effects of our larger dis- 
tance modulus (and assuming Mtot{UMi) = 2.3 x 10'' M©, 
after Mateo 1998) we obtain M/L ~ 7, just a factor ^ 4 
larger than the typical M/L of ordinary globular clusters. 

From the above discussion it appear that further ob- 
servational and theoretical efforts are still needed to fully 
understand the nature of UMi and Draco. It is quite pos- 
sible that we have some fundamental lesson to learn from 
these faint and unassuming stellar systems. 
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Bellazzini, Ferraro & Pancino 



Table 1 

Photometry and positions of stars in Ursa Minor. 
ID V €v I e/ Xjp^j Yjp^j 02000 "52000 Var. type Var. n. 



10018 


16.858 


0.005 


15.342 


0.005 


246.08 


365.12 


15" 


08™ 


27.2' 


67° 


10' 


07.7" 








10019 


17.018 


0.006 


15.628 


0.006 


399.38 


462.34 


15" 


08™ 


34.4' 


67° 


10' 


34.2" 








10020 


16.542 


0.004 


15.717 


0.005 


754.09 


801.17 


15" 


08™ 


51.2' 


67° 


12' 


07.0" 








10021 


16.579 


0.005 


15.620 


0.008 


1959.44 


817.81 


15" 


09™ 


48.4' 


67° 


12' 


09.3" 








10023 


16.928 


0.004 


15.679 


0.006 


838.95 


1488.32 


15" 


08™ 


55.4' 


67° 


15' 


15.9" 








10024 


16.458 


0.005 


15.569 


0.008 


873.74 


1918.16 


15" 


08™ 


57.2' 


67° 


17' 


14.2" 








10025 


17.714 


0.005 


16.489 


0.005 


548.78 


65.36 


15" 


08™ 


41.5' 


67° 


08' 


44.9" 








10026 


17.459 


0.006 


16.467 


0.009 


102.13 


160.05 


15" 


08™ 


20.4' 


67° 


09' 


11.5" 








10028 


18.090 


0.006 


15.754 


0.007 


1399.35 


279.70 


15'' 


09™ 


21.7= 


67° 


09' 


42.5" 









Note. — A sample of the photometrie eatalog. The complete catalog is available in ASCII format in the electronic edition 
of the paper. The variable stars are numbered after the nomenclature by N88 (colunm 11) and classified according the same 
authors (column 10: 1 = RRab; 2 = RRc; 3 = Anomalous Cepheids; 4 = unclassified variables). 



Table 2 

Photometry and positions of stars in Draco. 
ID V ev I £/ ^\px] ^[px] 0:2000 ^2000 Var. type Var. n. 



10069 


16.744 


0.004 


15.963 


0.003 


1295.30 


338.62 


17'' 20™ 


34.3'' 


57° 


51' 


56.4" 








10085 


17.510 


0.004 


16.742 


0.003 


1375.32 


999.17 


17" 20™ 


37.2' 


57° 


54' 


58.2" 








10086 


17.478 


0.003 


16.802 


0.003 


1145.29 


1051.62 


17" 20™ 


29.3' 


57° 


55' 


13.1" 








10090 


16.989 


0.004 


15.710 


0.002 


505.79 


1158.47 


17" 20™ 


07.2' 


57° 


55' 


43.5" 








10091 


17.313 


0.003 


15.744 


0.003 


1828.06 


1218.93 


17" 20™ 


52.9' 


57° 


55' 


57.8" 








10097 


17.112 


0.003 


15.816 


0.002 


463.32 


1303.23 


17" 20™ 


05.7' 


57° 


56' 


23.4" 








10100 


17.212 


0.004 


15.759 


0.002 


262.30 


1511.18 


17" 19™ 


58.8' 


57° 


57' 


20.9" 








10109 


17.193 


0.003 


16.272 


0.002 


1456.76 


1651.69 


17" 20™ 


40.2' 


57° 


57' 


57.7" 








10110 


16.948 


0.003 


15.716 


0.002 


270.91 


1652.71 


17" 19™ 


59.2' 


57° 


57' 


59.8" 









Note. — A sample of the photometric catalog. The complete catalog is available in ASCII format in the electronic edition 
of the paper. The variable stars are numbered after the nomenclature by Baade & Swope (1961) (column 11). In column 10 
only the known Anomalous Cepheid V141 is flagged with Var. type = 3. 



The Draco and UMi dwarf spheroidals 



Table 3 



Physical parameters for Ursa Minor and Draco. 





UMi 


Draco 


ref. 


E(B-V) 


0.03 


0.03 


M98 


Vrr 


19.86 ± 0.09 


20.28 ± 0.10 




VzAHB 


20.02 ± 0.09 


20.44 ± 0.10 




^bump 


19.40 ± 0.06 






(m - M)o 


19.41 ± 0.12 


19.84 ± 0.14 




rc 


17.9' 


7.7' 


P02, OOl 




77.9' 


40.1' 


P02, OOl 


e= (l-b/a) 


0.80 


0.29 


IH95 


My 


-10.3 ±0.4 


-9.0 ±0.3 


M98, P02, this work 


My(0)[mag/n"] 


25.5 ±0.5 


25.3 ±0.5 


M98 


< [Fe/H] > (ZW) 


-1.8 ±0.1 


-1.7±0.1 




[Fe/HUod (ZW) 


-1.9 ±0.1 


-1.6 ±0.1 




< [Fe/H] > (CG) 


-1.6±0.1 


-1.5±0.1 




[Fe/HUod (CG) 


-1.7±0.1 


-1.4±0.1 






0.10 


0.13 




0.12 


0.15 





Note. — A summary of relevant observational parameters for UMi and 
Dra. If no explicit reference is indicated, the estimated obtained in the 
present analysis is reported. The reported My has been obtained adopting 
our distance moduli. Lcgcnda of the previously undefined acronyms: P02 
= Palma et al. (2002); OOl = Odenkirchen et al. (2001); IH95 = Irwin & 
Hatzidimitriou (1995). The reported values of My have been obtained from 
the apparent total magnitudes listed by M98 (corrected according to the 
results by P02, in the case of UMi) and the distance moduli obtained in the 
present analysis. 



